Spectrum and thermodynamic currents in one-dimensional Josephson elements 



A. Krichevsky*, M. Schechter, Y. Imry, and Y. Levinson 
Department of Condensed Matter Physics, The Weizmann Institute of Science, Rehovot, 76100, Israel 
(Submitted January 28, 1999, Published in Phys. Rev. B 61, 3723 (2000)) 

The dc Josephson effect is considered from the thermodynamic point of view. Universal ther- 
modynamic equations, relating both bound and continuum contributions to the Josephson current 
with the normal electron scattering amplitudes are derived for the single mode case. To derive these 
equations we use and further develop the method of spatial separation between the superconduct- 
ing and normal parts of the junction. We also use this method to find the Andreev bound states 
in structures containing superconducting components. The general thermodynamic formulas are 
applied to the calculation of the current in various Josephson-type structures. In particular, the 
crucial role of the continuum contribution is demonstrated, even for short junctions (where it is 
usually neglected). We also find structures where the bound states supporting the giant currents 
are well separated; thus they can, hopefully, be populated nonuniformly and such current can be 
measured. 



I. INTRODUCTION 

Numerous theoretical investigations of Josephson junc- 
tions show that the equilibrium nondecaying current can 
be expressed in terms of a quasiparticle description based 
on Bogolubov-de Gennes (BdG) equations. There are 
two contributions to the current: one is supported by 
discrete states lying within the superconductor gap (so 
called "bound" current); and the other, carried by contin- 
uum of propagating modes outside the gap ( "continuum" 
current). Most authors find the bound contribution from 
the thermodynamic relation 
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where F is the free energy of the system and tp is the su- 
perconductor phase difference across the junction, while 
the continuum one is found from a Landauer type con- 
sideration which yields a current resulting from imbal- 
ance of left-going and right-going quasiparticle fluxes (it 
is non-zero in the superconductors). 

A more general method, allowing to obtain both dis- 
crete and continuum contributions from the thermody- 
namic approach was suggested by Beenakker Jlj . He de- 
rived the following expression for the free energy: 
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where "const." represents the (^-independent term, can- 
celing the divergence of the first term at high energies. 
The only significant approximation made in the deriva- 
tion of this formula is the steplike pair-potential shape: 
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where L is the junction length. The applicability of this 
model was discussed previously 0-0] ■ The important 



conditions for applicability of this model are: (i) Exis- 
tence of a high barrier on the boundary of the super- 
conducting and normal media, (ii) Low current, such as 
hep s /(mJQ) (where Jq is the current density and p s - 
the superconducting density) to be much greater than 
all characteristic lengths such as coherence length, junc- 
tion width, etc. The superconductor - semiconductor - 
superconductor and superconductor - insulator - super- 
conductor junctions usually belong to this category. 

The expression for the Josephson current follows di- 
rectly from Eqs. (|]) and (g): 
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where p is the continuum density of states jjj . The first 
term of Eq. (^) is the well known bound state current 
while the second term represents the continuum 
contribution, ft requires the knowledge of the depen- 
dence of the continuum density of states (DOS) on ip. 
An important tool, providing this information is Krein's 
theorem (see Appendix A), that connects the change 
in the DOS induced by any scatterer with the correspond- 
ing scattering matrix. 

Using this theorem, Beenakker jlj expressed the con- 
tinuum contribution via the scattering matrix of the nor- 
mal electrons in the multichannel case. We are going to 
specialize this approach to one-dimension (ID) case de- 
veloping the simplified expressions for this special case. 
Below we consider a general junction, consisting of an 
arbitrary one-dimensional (ID) nonsuperconducting con- 
striction (barrier) sandwiched between two superconduc- 
tors. The relevant scattering matrix has the dimensions 
4x4 and it can be expressed in terms of the energy- 
dependent normal electron S-matrix and the scattering 
amplitudes on the NS boundaries. 
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FIG. 1. (a) General S — barrier — S and (b) 
/ — impurity — I — S structures. (1) Geometric outline of 
the real problem. (2) Our model problem. (3) The energy 
diagram in the latter case. Thick line represents the potential 
V(x), thin line: the pair potential, dashed line: the elec- 
tro-chemical potential. 



II. SPATIAL SEPARATION OF THE 
SUPERCONDUCTORS FROM THE BARRIER 

In the next sections we express the bound spectrum 
and the Josephson current in a superconductor/general 
nonsuperconducting constriction/superconductor (SCS) 
junction via the scattering matrix of normal electrons. 
To do so we insert a fictitious ideal clean normal leads of 
length I between the superconductors and the constric- 
tion [Q [Fig. 0(a)]. The propagation of electrons in the 
"normal" medium is described by the model Hamiltonian 
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Below we consider the limit I — ► 0. Beenakker when 
inserting the fictitious normal metal, demanded that its 
size is smaller than all lengths in the system except the 
Fermi wavelength, since in this case, he claimed, the nor- 
mal metal insertion does not affect the results. This is 
true in the zero order in A/fi (Andreev approximation). 
Some of our calculations are done beyond the Andreev 
approximation, and in this case one must take the length 
of the normal metal to be shorter than the Fermi wave- 
length. This would appear meaningless from the physi- 
cal point of view. However, we adopt the formal point 
of view, and consider the "normal metal" as a domain 
in the BdG equation in which the external potential and 
the order parameter are both zero. We then can take this 
domain to be shorter than the Fermi wave length. On the 
other hand, the scattering formulation is still valid even 
when the normal-metal length is infinitesimally small be- 
cause in one dimension the outgoing wavefunction takes 



its asymptotic form cx e tkx at any distance from the bar- 
rier, however small it is, since there are no evanescent 
modes in the ID case. 

Below we demonstrate that our results for the energy 
spectrum, persistent currents (when exist), reflection am- 
plitudes etc. for arbitrary structures with fictitious nor- 
mal layers agree with the corresponding results for analo- 
gous structures without these layers obtained previously 
using different methods (|^,^| both under Andreev ap- 
proximation (Appendix B) and beyond it (Sec. III). This 
is to be expected, since the insertion of the fictitious 
normal metal corresponds to a non-singular perturbation 
and it does not influence the results in the limit I — > 0. 

Where the Andreev approximation is used, it is con- 
venient to set the Fermi momentum in the normal and 
superconducting leads to be equal. In this case one can 
neglect the normal reflection at the NS boundaries H . 



III. BOUND STATES IN THE INIS JUNCTION 

To illustrate the spatial separation method we con- 
sider the formation of the bound states in the insulating- 
superconducting structure containing an " impurity" em- 
bedded in the insulating region. By impurity we mean a 
quasi-bound state of energy Eq and width T that would 
be formed in analogous system where the superconduc- 
tor is replaced by a normal metal. Such a situation takes 
place, for instance, in the INIS junction [Fig. 0(b)]. An 
analogous problem has been considered by Laikhtman in 
terms of the tunneling Hamiltonian approximation [ jTo| . 
It was found that there exists a single bound state in this 
system, whose energy can be found from the equation: 



T 2 - E 2 = 
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In the limit r — * the energy E of this state tends to 
Eq if —A < Eq < A. On the other hand, using more ex- 
act wavefunction matching method it was shown H that 
there exists a bound state of energy close to A at an IS 
boundary without impurity. 

We consider the following problem: An impurity with 
energy Eq, — A < Eq < A located at a distance L from 
the superconductor. Using the spatial separation method 
we find two energy levels: one close to Eq and one close 
to A for the case of large L (small T). The level close to 
A can not be obtained using the tunneling Hamiltonian 
method. 

In order to find the bound states in this structure we 
insert a fictitious infinitesimally short normal layer be- 
tween the superconductor (S) and the rest of the struc- 
ture (R). Consider the motion of the quasiparticle in the 
normal layer. The quasiparticle can be either normal or 
Andreev reflected at the NS boundary and only normal 
reflected at the RN boundary (see Fig. ||). 
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FIG. 2. Formation of a bound state on RS boundary. Solid 
and dashed arrows represent the electron and the hole prop- 



agation correspondingly. Amplitudes 
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the Andreev and normal reflection amplitudes for electrons 
(holes) on NS boundary; r' e , h \: the normal reflection on RS 
boundary (prime indicates that the reflected particle moved 
initially to the left), numbers enumerate the scattering chan- 
nels. 



From the uniqueness of the wave function one obtains 
the eigenenergy equation: 
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In the case that R contains an impurity we assume that 
the reflection amplitude of electrons (holes) at energies 
close to the bare impurity level take the form 



r e (E) =r (E) 
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that contains both a Breit-Wigner-type resonance and 
a slow energy dependence ro(E) as in Eq. (|^) below. 
Note, that \r e \ = 1 since there are no propagating modes 
within the insulator. Substituting these reflection am- 
plitudes to Eq. (Q), neglecting the energy dependence of 
ro and using the Andreev approximation which provides 
b e = b% = 0; a e = ah = vq/uq [where uq and vq are 
the BCS coherence factors: u§ = (1 + y/l - A 2 /E 2 )/2, 
v$ = (1 - y/l - A 2 /E 2 )/2], one obtains Eq. (§) which 
has only one positive energy solution. Taking into ac- 
count the energy dependent ro (E) is crucial for obtaining 
the second level. 

The numerical solution of Eq. (0) for the case /i = 
100A, V = 130A, En = 0.5A [see Fig. 0(6)jand the re- 
flection amplitudes defined in Eqs. (H) and (£)]) (below) is 
shown in Fig. |^. The exact calculations for INIS struc- 
ture (rather than Breit-Wigner-type approximation (j^) 
give analogous results. The Breit-Wigner-type approx- 
imation is used here just for a clear definition of the 
resonance width T and for direct comparison with the 
tunneling Hamiltonian method. 
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FIG. 3. Typical behavior of discrete levels as a function 
of the resonance width for the Breit-Wigner type reflection 
amplitudes. The higher level disappears in the continuum at 
very small V, so the upper left part of the graph is shown in 
the inset. 

For small T there are two bound states: the upper level 
lies close to the gap edge (its wave function is localized 
around the interface — mostly in the superconductor ) 
while the lower state is localized at the impurity and has 
an energy close to Eq. As the coupling between the im- 
purity and superconductor increases, so does the overlap 
between the bound-state wave functions causing the re- 
pulsion of the levels. The upper level soon disappears 
in the continuum. The energy of the lower level first 
decreases, but then starts to increase and tends to A 
asymptotically. 

For the upper (surface) level one can rederive, us- 
ing our formalism the result obtained by Wendin and 
Shumeiko || for IS boundary (they used the matching 
of the decaying wave functions on the boundary) by tak- 
ing the limit r — * 0. In this case, according to Eq.(||) 
r(E) — > r (E) as for pure IS boundary with E — E > T. 
We model the insulator by a potential step of height V, 
then the reflection amplitudes (see Fig. ||) are 
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where 



y/2m(n±E)/h; 
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("+" corresponds to electrons, "— " to holes). Substi- 
tuting these values of r e ,r^ and the Andreev reflection 
amplitudes calculated in the next section to Eq. (f^) and 
then expanding up to first order in A//i one finds the 
energy of the interface level: 
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Note, that expansion up to zero order in A/// leads to a 
wrong result {A - E = A 3 /[2(i(V - fj,)]}. This indicates 
that the Andreev approximation is to be used with care. 

Taking into account the energy dependence of ro in 
Eq. (||) is crucial in order to obtain the edge bound state. 
Neglecting the energy dependence of ro results in obtain- 
ing just one (impurity) level, exactly coinciding with the 
level obtained using the tunneling Hamiltonian method. 

The existence of such bound states helps to understand 
the discrete spectrum of Josephson elements. The SIS 
junction can be considered as two coupled IS boundaries 
H and the bound spectrum consists either of a pair of 
levels lying close to the gap edge (the levels determined 
by Eq. (|1J) split by coupling analogously to the levels in 
a double- well problem |ll[]) or of a single level when the 
coupling is strong enough for pushing the upper level to 
the continuum. Analogously, if the junction barrier has 
a transmission resonance within the gap (for instance, a 
SINIS junction or an element including the impurity), 
one can consider it as an "IS boundary with a resonance" 
coupled to the "pure IS boundary" and expect at most 
three discrete levels: one close to the resonance and two 
levels close to the gap edge if the coupling is weak; two or 
one level for stronger coupling. For N resonances within 
the gap there are possible N + 2,N + l,N . . . levels, de- 
pending on coupling strength. 

All these levels can be slightly modified by changing 
the phase difference between the superconductors, giv- 
ing rise to the bound state contribution to the Josephson 
current. Although the edge states often lie almost indis- 
tinguishably close to A, their contribution to the current 
can be major M. 
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FIG. 4. Simple illustration of scattering processes. Solid 
arrows represent electronlike waves, dashed arrows: holelike 
ones. Numbers enumerate the scattering channels. 

Second step: \ini) scatters on the barrier and converts 
to \0ut2}' 
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Here we used the symmetry S n providing t' = t and also 
the time-reversal symmetry of Schrodinger equation, as- 



suring 



r h {E) = r* e (-E); t h (E) = t* e (-E). 



(16) 



IV. S MATRIX 



In this section we obtain the scattering matrix of a 
general nonsuperconducting constriction sandwiched be- 
tween two superconductors. The knowledge of the scat- 
tering matrix and Krein's theorem, which relates the 
change in the continuum density of states to the S ma- 
trix, allow us to calculate the Josephson current of the 
continuum (next section). A simple way to obtain the S 
matrix is to reduce the whole scattering process to the 
scattering of normal electrons on the constriction and to 
the processes on pure SN boundaries. 

To simplify the discussion we divide the scattering pro- 
cess into several steps (Fig. ^): 

First step: An incoming wave \in) hits the SN bound- 
aries from the outside. A part of it penetrates into the 
normal leads forming the new "incoming wave" \in±) and 
and the rest is reflected as \out\): 



Third step: \0ut2) reaches the NS boundaries from in- 
side. Partially it penetrates into the superconducting 
electrodes forming the outgoing wave \0ut3) and partially 
it is reflected back to the normal region as \iri3): 



\in 3 ) = S A \out 2 ); \out 3 ) = c NS \out 2 ), 



(17) 



where the matrices Sa,cns are the Andreev reflection 
and transmission matrices analogous to a s ,csN- 

Fourth step: Continue the steps 2 and 3 up to infinity. 



The total output can be found from Eqs. (|13|),(14) and 
p7|) by summing up the geometric series: 

I out) = \outi) + \0ut3) + \0ut5) + . . . 

= [a a + c NS (I - SnSa^Sncsn]]™) = S\in) (18) 
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FIG. 5. The illustration of the scattering processes. Letters 
N, S correspond to the normal and superconducting media, 
the order parameter on both sides is shown. Solid arrows 
show the electronlike wave propagation, dashed ones show 
holelike ones. The index e means that the incoming parti- 
cle is electronlike, subscript NS means that it comes from 
the normal part (SN: from the superconductor; prime indi- 
cates that the particle comes from the right. The scattering 
processes for holes are completely analogous, one has just to 
replace the index e by h and solid arrows by dashed ones (and 
vice versa). 



The matrix Sn is given by Eq. (15) and the matri- 
ces (1 s ,csn,cns, Sa are constructed from Andreev and 
normal reflection and transmission amplitudes on SN 
boundaries which we obtain below. 

In their pioneering work Blonder, Tinkham, and Klap- 
considered the following scattering process (see 
5): the electron, propagating from left to right in the 
normal metal is scattered by the NS boundary. It can be 
reflected back to the normal medium or transmitted to 
the superconductor either as an electron or as a hole. We 
follow their treatment, but use a different normalization 
of the wave functions to assure unitarity of the scattering 
matrix. The incoming, reflected and transmitted waves 
for NS boundary are given by: 
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where 



fc* = \J2m(Li±v / E 2 - A 2 )/h. 



(20) 



Our normalization differs from the one used by Blon- 
der, Tinkham, and Klapwijk by prespinor factors 



The amplitudes a, b, c, d have been found by matching 
the wave functions and their derivatives on the bound- 
ary. We define the scattering amplitudes analogously 
for the processes on both NS boundaries of the model 
Josephson junction (Fig. |J) when the phase difference 
ip is maintained between the superconductors, consider- 
ing electronlike and holelike projectiles coming from the 
normal medium to the superconductor (corresponding 
amplitudes are indicated by NS below) and vice versa 
(index SN) — see Fig. ||. In contrast with the work of 
Blonder, Tinkham, and Klapwijk || we do not initially 
neglect the deviations of projectiles momenta from Uf- 

In terms of these amplitudes the scattering matrices 
have the following form: 
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(22) 



The scattering amplitudes are again found from the so- 
lution of the matching problem on the boundaries. Here 
we list them in explicit form: 



2Vk-k+{q + + q~ )u v a /D, (23) 
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-v 2 (q+ -k-)(q~ +k+)]/D, (24) 

c esN = 2^q+(q- + k-)u ^u 2 -v 2 e-^^/D, (25) 

d esN - 2^q~(k- - q+)v^u 2 -v 2 e^/ A /D, (26) 

a eNS = 2^q+(k+ + k-^voe-^^/D, (27) 
be NS = [u 2 (q + - k+)(q- + AT) 

~v 2 (q- ~k+)(q+ + k~)]/D, (28) 

c eNS = 2^/fc+^( ( ^ + k-)u^u 2 -v 2 e-^/ A /D, (29) 

d eNS = 2^k~q^(k+ - q-)v Q ^Jul-vle-^' A /D, (30) 



where 
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D = ul(q+ + k+)(q- + AT) - v 2 (q+ - k~)(q- - fc+). 



(31) 



Note, that the nonprimcd amplitudes correspond to 
the processes where the projectile comes from the left, so 
the above values a esN , b esN . . . are calculated on the left 
SN boundary (the phase of the order parameter is equal 

are on the 



to —ip/2), while the amplitudes a, 
right boundary (the phase of the order parameter ip/2). 

To obtain the expressions for hole scattering ampli- 
tudes (such as cih SN ) one has to replace in the above 
formulas , q ± by k^ , q^ correspondingly and also to 
replace ip by —<p. All primed amplitudes are equal to the 
corresponding nonprimed ones taken at —(p. Inserting 
these matrices and S n ([IB]) into Eq. ([l9|) we obtain the 
S-matrix for a general SCS junction. 



V. JOSEPHSON CURRENTS 

Knowing the explicit form of the scattering matrix one 
can find the ^-dependence of the continuum density of 
states using Krein's theorem: 
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The idea of the proof is given in Appendix A. 

Now we have all the required tools to find the Joseph- 
son current. The continuum contribution can be ex- 
pressed in terms of the scattering matrix using Eqs. (^) 
and (H) |@: 

E 

dE tanh — — — In det S. (33) 
dip 

This result can be simplified by using the Andreev ap- 
proximation (that is, neglecting the deviations of k ,q 
from fcp): 
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where 



D = \t(E)t(-E)\; = M gi(£) - argi(--E); ( 35 ) 
R=\r(E)r(-E)\; e = argr(£) - argr (-E) - f3; (36) 
cos 7 = R cos e + D cos ip (37) 

(we follow the definitions of Wendin and Shumeiko p|). 

Equation (|34]), which is one of the central results of our 
work, embraces a wide class of Josephson junctions. In 
particular, Eq. ( |34"| ) describes arbitrarily long junctions, 
a problem often avoided. Equation ( |34"| ) is convenient for 
qualitative estimations of the continuum contribution to 



the current. For instance, consider a clean SNS junc- 
tion. In this case (3 = EL/ (A£q) where £o = ^/{Akp) is 
the superconductor coherence length (for more details see 
Appendix B). For a short junction (L <C Co) Ic is small 
and it is roughly proportional to the junction's length. 
In the opposite case of a very long junction (L 3> £ ) the 
continuum contribution is small again, but for a different 
reason: the integrand is now a product of the oscillat- 
ing {sin[EL/(A£o)]} and decaying (« l/E 2 ) functions of 
energy and the integration over energy results in strong 
cancellation. However, in some junctions one can expect 
a significant continuum current: (i)If /?, 7 are small (as 
in the case of some SIS junctions [§) or, more general, 
if ± 7 is close to irn — the denominator in Eq. (H) is 
small at E > A. (ii) If the junction transmission D is 
large or if it has a sharp maximum at E > A (as in the 
case of resonant structures like SINIS, SININIS, etc.) 
and thus the numerator can also be large. 

If the scattering amplitudes of the constriction depend 
weakly on energy, D,R are approximately equal to the 
transmission and reflection coefficients of the barrier cor- 
respondingly and the angles /3, e are small. In this case, 
the integral in Eq. (Q) can easily be done at zero tem- 
perature: 
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The spatial separation method can also be used to red- 
erive the known results for the bound state energies and 
current. The eigenenergies E n correspond to the singu- 
lar points of energy-dependent scattering matrix. From 
equation ( p"9| ) we obtain B|: 



det[7 - S N (E n )S A (E n )] = 0. 



(39) 



Using the definitions of S N , S A [Eqs. @ and @], An- 
dreev approximation, properties (|16| ) of Sn and the pa- 
rameters -R, D, (3, e, 7 [Eqs. (35)-(p7pj one obtains the fol- 
lowing equations for the eigenenergies of the bound states 
in two equivalent forms: 
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+ D 2 cos ip, 



2 7 ^y / E 2 A 2 - E 4 sin f3 = R cos e + D cos ip. 
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Equation © can be solved with respect to E by squar- 
ing and some other algebraic manipulations. The result 
is 



E 2 



A 2 cos 2 if sin(/3 ± 7) > 0. 



(42) 



The condition sin(/3±7) > in ( [42] ) is essential, because 
squaring can produce a redundant solution [tt3| . 
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FIG. 6. Exact numerical calculations for the symmetri- 
cal SINIS junction, T ~ 0.04A. Dashed line represents the 
low-energy state current, dash-dotted shows the edge state 
current, "+": the total bound state contribution, "x": the 
continuum contribution, solid line: the total Josephson cur- 
rent. 

Although the formulas (|l|) and @) look like ex- 
plicit equations for the energy, there might be hidden 
energy dependence of the scattering amplitudes and con- 
sequently of the angles /3, 7, e. Thus equations (41) and 
( [42l ) are to be solved self-consistently. 

Knowing the scattering properties of the barrier one 
can calculate the (^-dependent current using Eq. ([34]) and 
one of Eqs. (pp|) , (|4l|) , and ([42]). Below we consider several 
applications of these formulas. 



VI. SINIS JUNCTIONS 



For the structures containing two IS boundaries (like 
SIS, SINIS etc.) one can safely use the steplike pair 
potential approximation, |^,[l4j] and consequently Eqs. 
(H),@),(il]), and (H||) for quantitative calculations. 
The SINIS junction is the simplest example of the res- 
onant structure. We solve it analytically using the Brcit- 
Wigner-like formulas for transmission and reflection am- 
plitudes. More precise calculations, taking the exact am- 
plitudes t, r instead of the Breit-Wigner approximation 
were performed numerically. 

Consider the symmetrical SINIS junction. Denote 
the TV-part length by L and the single-barrier transmis- 
sion coefficient by D. The Breit-Wigner near-resonance 
transmission and reflection amplitudes are 



* = 



r 



E — En 



E -E + iT 



E- E + iT 



, (43) 



where Eq is the bare resonance level, T = HDvf/L and 
(fit , fr are some weak (on the scale of T) energy depen- 
dent phases which cancel in further calculations. 



Currents, eT(E=0)A/fi 




FIG. 7. Symmetrical SINIS. The resonance width and all 
line types are the same as at the previous figure, but the res- 
onance itself is not close to the middle of the superconductor 
gap. The total current is reduced by an order of magnitude 
compared to the previous case. 

Substituting the Breit-Wigner amplitudes from 
Eq. (||) to Eq. @) we obtain 



E 2 {El + T 2 + A 2 ) + 2REVA 2 - E 2 



E 4 



£ 2 A 2 



A 2 r 2 cos 2 ^. 

2 



(44) 



This equation has only one positive energy solution and 
thus it misses the edge states as a result of inadequate 
accuracy of the Breit-Wigner approximation. Equa- 



tion (44) is accurate enough only in the neighborhood 
of the resonance, but not on the scale of A. In order to 



obtain the edge states in Sec. |I| we have introduced the 
slow energy dependence of rn. Now assume that the res- 
onance is sharp, r <C A. One can see [for example from 
the graphic solution of Eq. ([|I])] that when the bare reso- 
nance level lies within the gap, the bound state is close to 
|-Eq|- In particular, for £ <C A equation ( ]44| ) simplifies 
1 to 



E = 



El + T 2 cos 



(45) 



From Eq. ( J45| ) it follows that the amplitude of the <p de- 
pendence of such levels is of order of Y 2 /En when En T 
and of order of T for Eq <g; T. If the system allows to 
tune the bare resonance position one can observe the en- 
hancement of the current whenever Eq crosses the Fermi 
level |5Ul5|l. In the latter case E = T cos tp/ '2. 



In Sec. Ill we saw that for structures containing a single 
superconductor there exist bound states on the resonance 
as well as on pure IS boundary. For SINIS junction one 
could expect an existence of three bound states: one close 
to the l-Erj I and two others close to A. Exact numerical 
calculations not using Breit-Wigner approximation but 
exact amplitudes t, r for SINIS junction confirm this 
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expectation. The contribution of the "edge" states is, 
as in the SIS case, proportional to off-resonant trans- 
mission, that is to r 2 . Therefore, these levels can not be 
found using the Breit-Wigner approximation. Sometimes 
one of the "interface" levels (or both of them) is pushed 
to the continuum. As for the continuum contribution, it 
is proportional to the junction transparency (|34|) which 
is of the order of T 2 . 

The numerical solution for the case of a very low lying 
resonance (Eq -C T) is shown in Fig. |(| It is in a good 
agreement with the Breit-Wigner model discussed above. 
However, neither the edge state nor the continuum con- 
tribution can be neglected when Eq > T. In the example 
shown in Fig. [?] one edge level is pushed to the continuum 
(as in the previous case), but the other one contributes 
and together with the continuum reduces the total cur- 
rent by an order of magnitude. A similar reduction was 
observed by H.Takayanagi (private communication). 



VII. SININIS JUNCTIONS 

The SININIS junction represents two coupled nor- 
mal layers. For a short symmetrical junction without res- 
onances one can find "giant currents" similar to the ones 
m SNINS @. These are the contributions proportional 
to the square root of the junction transmission rather 
than to the transmission itself. However, the current- 
carrying levels are close to each other and thus they are 
almost equally populated in equilibrium and the "giant" 
currents strongly compensate each other. This results in 
a total current which is approximately equal to 



standard 



eAD(E = 0) 



(46) 



a well-known Josephson result. It is interesting to find 
a structure where the current exceeds this "standard" 
value. 



Junction transmission 



T(H) 



Currents, eT(E=0)\/fi 




0.5 10 0.5 1 

FIG. 8. The "twin resonance" structure. The bound and 
the continuum currents have the same sign. 



In this section we present the numerical results for two 
structures with resonances. 

The first one is the short symmetrical SININIS with 
twin resonances. We have chosen the constriction with 
one of twin peeks placed within the continuum (Fig. |^). 
Interesting, that the continuum current exceeds the total 
bound state current and that both discrete and contin- 
uum contributions have the same sign (they often have 
opposite signs). 

The second example (Fig. ||) is the long (about 3 
£o) nonsymmetrical junction. One can see again the 
strong compensation of individual levels, but the total 
current exceeds the standard value eAD(Q)/h by about 
a factor of 2. More important, the individual level cur- 
rents exceed the J^ tandard hundreds of times and the 
levels, supporting these currents are split significantly 
(E 1 w 0.75A,£ 2 ~ 0.89A,i;3 rs 0.95A). Thus it may 
be possible to populate the bound states nonuniformly 
(for example by resonant electromagnetic pumping or by 
coupling to the additional electrode ||) and to enhance 
the Josephson current by order(s) of magnitude. 



Junction transmission 



T(E)0-08 
0.06 
0.04 



Currents, eT(E=0)A/?j 



, , I 

-2-1 1 2 

Total current, eT(E=0)A/?j 

Of 
-0.2 
-0.4 
-0.6 



^bound 

+ + 

I, 



FIG. 9. Currents in a long (~ 3£o) nonsymmetrical 
SININIS junction. In this case there are three bound states 
and, correspondingly three bound currents (II, 12, 73). Al- 
though each one of them exceeds the "standard" value for 
Josephson current many times, the total current is just two 
times greater than that. 



VIII. CONCLUSIONS 



In the present work the dc Josephson effect was in- 
vestigated and explicit formulas for Josephson current 
were obtained in ID case. We used the approach sug- 
gested by I.O.Kulik |Q based on the solution of Bogol- 
ubov - de Gennes equations. This technique treats the 
contributions of discrete and continuum energy states 
separately. Unlike the usual way |^,^,|l^], we have not 
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used the Landauer-type consideration to find the con- 
tinuum current, but Beenakker's approach allowing to 
derive both discrete and continuum contributions from 
the most general thermodynamic relations. We have also 
specified the Beenakker's idea of spatial separation of the 
superconductors from the barrier for ID case and devel- 
oped it for infinitesimally thin separating layers (they do 
not have to be long compared to the Fermi wavelength 
Xf)- Such insertion of fictitious normal layers should be 
treated as a mathematical trick only, reducing quite com- 
plicated Josephson problem to Andreev reflection PJlTj 
and relatively simple scattering problem of normal elec- 
trons on the barrier. 

Our results are applicable to a wide class of ID ele- 
ments. The only approximations we used are: (i) Andreev 
approximation, A -C /i. (ii) Existence of a large barrier 
on the boundary of the superconductor, (iii) Low cur- 
rent, heps/ (mJo) (where Jq is the current density and 
p s is the superconducting density) to be much greater 
than all characteristic lengths such as coherence length, 
junction width, etc. The last two conditions are required 
to justify the steplike pair potential approximation (^). 

Applying our formulas to different junctions we have 
appreciated the crucial role of the continuum contribu- 
tion even in cases where this was not expected. The con- 
tinuum current is often neglected in short (compared to 
the coherence length) junctions 0j^], but in some short 
structures we found that the continuum current can be 
of order of the bound-state current or even exceed it. We 
have found also an unusual enhancement (rather than 
decrease) of the total current by the continuum contri- 
bution. 

It was found in the paper by Wendin and Shumeiko 
that total Josephson current might result from almost 
complete cancellation of huge bound state currents flow- 
ing in the opposite directions. The states supporting 
these giant currents lie very close to each other, so it is 
hard to populate them differently and the giant currents 
are almost canceled. We found some structures which we 
believe can be built experimentally, where the level sep- 
aration is of order of A and the individual currents sup- 
ported by these states exceed the usual Josephson value 
eAD(E = 0)/h tens or hundreds of times. We still could 
not avoid the strong cancellation of individual currents 
and the equilibrium current is of order of its standard 
value (see above), and in our best structures it is en- 
hanced by a factor of 2-3. We expect that the current 
in such junctions can be much enhanced by appropriate 
population of individual bound states (for example, by 
microwave pumping or by coupling to another electrode) . 
However, just a large energy separation of Andreev states 
does not guarantee the lack of cancellation of their con- 
tributions. 

We benefited from stimulating discussions with B. 
Laikhtman and V. Shumeiko. We also thank E. 
Balkovsky for his kind help. This work was supported 



by the Israel Academy of Science and by the German- 
Israeli Foundation (GIF). 



APPENDIX A: KREIN'S THEOREM 

In this appendix we find the relationship between the 
continuum density of states to the scattering matrix. We 
use the definition of the density of states for the contin- 



p = Im Tr G 

7T 



(Al) 



where G is the retarded Green's function G = (E — H + 
ie)- 1 , H is considered to be a full Hamiltonian of the sys- 
tem, including the scatterer contribution: H = Hq + V, 
and Hq is an unperturbed one. In terms of these variables 
the scattering operator takes the form 



S — = Gr\ GGGr\ 



(A2) 



Here ft± = [E - H ± ie] ~ 1 [E - H ± ie] is the Moller wave 
operator. Krein's theorem JjJ claims that for any two 
linear operators H± , Hi (for example for the free Particle 
Hamiltonian and the perturbed Hamiltonian) holds: 



Tr[{H x - EI)- 1 - (H 2 - EI)- 1 ] 
d 



dE 



lndet[(#i - EI)(H 2 - EI)- 1 }. 



(A3) 



The idea of the proof is simple: in the basis of eigen- 
functions the operator takes the diagonal form, so 



d_ 

oe 



lndet(-ffi — EI) 



8 



( E\ — E 



dE 

= -E 



In det 



Eo-E 



E n - E 



-Tr(i?i - Eiy 



(A4) 



The theorem stated follows from Eq. (A4) and the anal- 
ogous relation for H^. Note that it is not required to 
diagonalize both Hamiltonians simultaneously. From the 
definitions of the DOS (Al) and the definition of the 



scattering matrix (A2) one obtains 



1 9 , , „ 
2m oE 



(A5) 



where po is the free particle DOS and p is the DOS of 
the scattering problem. We also used the unitarity of the 
scattering matrix. 

The theorem is valid in both the normal and the su- 
perconducting cases, in the latter case one has to treat 
the corresponding Hamiltonians as BCS or Bogolubov 
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Hamiltonians and to use the superconductor - to - su- 
perconductor scattering matrix. Noting that po is if in- 
dependent (up to small mesoscopic corrections) one can 
calculate the continuum contribution to the current using 
Eqs. (§) and (|A5|). 



matching problem (8x8 linear system). Its advantage be- 
comes even more pronounced in more complicated prob- 
lems with larger number of boundaries. In the limit I — > 
our results tend to the ones obtained by Blonder, Tin- 
kham, and Klapwijk Q for reflection from NS boundary 
with a barrier [Eg). 



APPENDIX B: APPLICATIONS OF THE 
SPATIAL SEPARATION METHOD 



2. SNS junction. 



1. Andreev reflection 

In this appendix we rederive the formulas for Andreev 
reflection from NS boundary with a barrier using the 
spatial separation method. Consider the reflection of a 
quasiparticle from a point "impurity" modeled by a 5- 
function barrier, separated by a distance I from an ideal 
NS boundary. More precisely, we use the Bogolubov- 
de Gennes [|l4| Hamiltonian with V(x) = +VoS(x) and 
A(x) — A 9(x + I) where 9 is the Heaviside step func- 
tion. Taking into account multiple reflections from the 
barrier and Andreev reflections from ideal (barrier free) 
NS boundary H|TjJ we obtain for incoming electron-like 
particle 



1 - r' f: r h a e a h e 2i( -i+-<i-) 1 

r' h tla e a h e 2i ( q+ - q ~V 
1 - r' e r' h a e a h e 2i( -i + ^^^ 



(Bl) 
(B2) 



where a, b are the total Andreev and normal reflection 
amplitudes, r e ,t e ,rh,th are the barrier reflection and 
transmission amplitudes for electron and hole, a e , ah are 
analogous Andreev reflection amplitudes, prime corre- 
sponds to the left-going particle. Restricting our con- 
sideration to energies of order of Ao and neglecting all 
Ao//i terms (Andreev approximation) we have 



a e — a/i — — j 

M 



r e . = r = 



-iZ 
1 + iZ ' 



t, = t = 



1 



1 + iZ 
(B3) 



where Z = 2mVo/(h 2 kF)- Corresponding normal re- 
flection and transmission amplitudes for holes are just 
the complex conjugated ones for electrons. Substituting 
these quantities to Eq. (B2) we obtain 



u v e 



b=- 



u 2 + {u 2 -v e 2 ^ + -i-) l )Z 2 ' 
VQe ~i{i + -i-)l)(iZ + Z 2 ) 



(U 2 - ^ P 2 ^ 



ug + (u 2 - voe^-i-^Z 2 



(B4) 
(B5) 



This intuitively clear two-step method allows to get the 
final result more easily than by a direct solution of the 



As an illustration of the previous results, consider the 
well-investigated example of SNS contact (here N indi- 
cates a normal metal, S a superconductor). 

The SNS constriction includes no barrier, it consists 
of a piece of a clean normal metal of length L sandwiched 
between the superconductors. For this structure 



e -iq l . f) — i 

(B6) 



r(E) = 0; t(E) = e l " +L ; t*(-E) 



Substituting these expressions to Eq. ( |40[ ) and using the 
relation vq/uq — e ~ tm ' ccosE / A we obtain the Kulik's re- 
sult for bound states [0: 



2 arccos ■ 



E n 
A 



(in -q n ) L±t P= 27rn - 



(B7) 



The momenta q + ,q are very close to fc^- Expanding 
q + , q~ as a function of energy around Uf we obtain 



k F ± 



E 



2A£ 



(B8) 



where £o = / i /(A)(l/fcF) is the superconductor coher- 
ence length. Using Eqs. (|34|), (B6), and flB§| ) one can 
find the continuum contribution: 



T A 2 

I c = —A smcp 

7m 

EVE 2 -A 2 sin (EL/A£ ) tanh (E/2k B T) dE 

x (E2_ A 2 COS 2 {[(EL/Ai )- V ]/2})(E 2 -A 2 cos 2 {{(EL/AZ )+ v ]/2}) • 

An analogous equation was obtained by Bagwell (l(J . Re- 
lations (B7) and (B9) should not be treated too seriously 
in quantitative aspect. The problem is that the step-like 
pair potential hypothesis fails; conversely, A changes on 
scale £o 0. 

* Present address: University of Alberta, Edmonton, 
AB, Canada T6G 2J1. 
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